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Detecting  Spectrum  Opportunities  in  Poisson  Primary  Networks 

Wei  Ren  and  Qing  Zhao 


I.  Introduction 

Opportunistic  spectrum  access  (OSA),  also  referred  to  as  spectrum  overlay,  is  one  of  the 
several  approaches  envisioned  for  dynamic  spectrum  management  [1],  The  basic  idea  of  OSA 
is  to  allow  secondary  users  to  exploit  temporarily  and  locally  unused  channels  without  causing 
unacceptable  interference  to  primary  users. 

One  of  the  central  issues  in  OSA  is  spectrum  opportunity  detection  through  sensing.  Before 
transmitting  over  a  particular  channel,  a  secondary  user  needs  to  decide  whether  this  channel  is 
an  opportunity.  This  is  the  so-called  “Listen-before-Talk”  (LBT).  In  this  report,  to  quantitatively 
characterize  the  performance  of  LBT,  we  consider  a  slotted  network  setup  where  primary  users 
are  distributed  according  to  a  spatial  Poisson  process  (see  Sec.  II).  We  first  analyze  the  perfor¬ 
mance  of  LBT  at  both  physical  and  MAC  layers,  and  then  we  investigate  the  translation  from 
physical  layer  opportunity  detection  performance  to  MAC  layer  performance,  which  is  crucial 
in  the  design  of  the  opportunity  detector.  In  particular,  we  demonstrate  the  complex  dependency 
of  the  relationship  between  PHY  and  MAC  on  the  application  type  (guaranteed  delivery  vs. 
best-effort  delivery)  and  the  use  of  handshaking  signaling  such  as  RTS/CTS  at  the  MAC  layer. 

Throughout  the  report,  we  use  capital  letters  for  parameters  of  primary  users  and  lower-cased 
letters  for  secondary  users. 


II.  Network  Model 

Consider  a  decentralized  primary  network  with  slotted  transmission  structure.  Assume  that 
users  are  distributed  according  to  a  two-dimensional  homogeneous  Poisson  process  with  density 
A.  At  the  beginning  of  each  slot,  each  primary  user  has  a  probability  p  to  transmit  data  to  a 
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receiver  that  is  uniformly  distributed  within  its  transmission  range  Rp  (see  Fig.  1).  Based  on 
the  Thinning  Theorem  and  the  Displacement  Theorem  for  marked  Poisson  processes  [2],  both 
primary  transmitters  and  receivers  form  a  two-dimensional  homogeneous  Poisson  process  with 
density  pX.  Note  that  these  two  Poisson  processes  are  not  independent. 


Fig.  1.  A  Poisson  distributed  primary  network. 


III.  Spectrum  Opportunities:  Definition  and  Implications 

A  rigorous  study  of  spectrum  opportunity  detection  must  start  from  a  clear  definition  of 
spectrum  opportunity.  In  fact,  the  concept  of  spectrum  opportunity  is  more  involved  than  it  at 
first  may  appear  [3].  As  illustrated  in  Fig.  2,  a  channel  is  an  opportunity  to  a  pair  of  secondary 
users  A  and  B  if  they  can  communicate  successfully  while  limiting  the  interference  to  primary 
receivers.  More  specifically,  B  can  not  be  affected  by  primary  transmitters,  and  A  can  not 
interfere  with  any  primary  receivers.  In  other  words,  there  is  no  primary  receiver  located  within 
distance  ry  to  A,  and  no  primary  transmitter  within  distance  Rj  to  B,  where  77  is  the  interference 
range  of  secondary  users  and  is  monotonically  increasing  in  the  transmission  power  ptx,  Rj  is  the 
interference  range  of  primary  users.  Here  we  have  adopted  the  disk  model  for  signal  propagation 
and  interference.  This  definition,  however,  applies  to  general  cases  [5]. 

The  first  implication  of  the  spectrum  opportunity  is  that  to  detect  the  opportunity,  A  needs 
to  detect  the  presence  of  nearby  primary  receivers,  and  B  needs  to  detect  the  presence  of 
nearby  primary  transmitters.  However,  without  assuming  cooperation  from  primary  users,  primary 
receivers  are  much  more  difficult  to  detect  than  primary  transmitters.  One  effective  alternative 
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Fig.  2.  Definition  of  spectrum  opportunity. 


is  ”listen-before-talk”  which  transforms  the  problem  of  detecting  primary  receivers  to  detecting 
primary  transmitters.  As  shown  in  Fig.  3,  A  infers  the  presence  of  primary  receivers  within 
distance  ry  from  the  presence  of  primary  transmitters  within  distance  ro,  where  ro  denotes 
the  detection  range  and  can  be  adjusted  by  changing,  for  example,  the  threshold  of  an  energy 
detector.  Let  Rp  denote  the  transmission  range  of  primary  users,  then  the  most  conservative 
detection  range  is  Rp  +  rj,  which  may  result  in  overlooked  opportunities.  As  illustrated  in  Fig  3, 
if  the  detection  range  is  Rp  +  rj,  the  transmission  activities  of  primary  transmitter  X  will  prevent 
A  from  accessing  the  spectrum  opportunity  even  though  the  intended  receiver  of  X  is  outside 
the  interference  range  ry  of  A.  However,  if  we  let  m  <  Rp  +  rj,  e.g.,  the  dashed  circle  in  Fig.  2, 
A  can  not  detect  the  transmission  activity  of  Y,  and  it  will  interfere  with  the  intended  receiver  of 
Y .  So  even  if  A  uses  LBT  with  ’’perfect  ears”,  spectrum  opportunity  detection  can  not  achieve 
perfect  performance  due  to  the  fundamental  deficiency  of  LBT. 

The  second  implication  is  that  spectrum  opportunity  is  asymmetric,  that  is,  a  channel  that 
is  an  opportunity  when  A  is  the  transmitter  and  B  the  receiver  may  not  be  an  opportunity 
when  B  is  the  transmitter  and  A  the  receiver.  So  if  B  needs  to  send  back  the  acknowledgement 
to  A  after  B  receives  data,  the  acknowledgement  from  B  may  not  be  received  by  A.  This 
asymmetry  leads  to  a  complex  relationship  between  the  opportunity  detection  performance  at 
the  physical  layer  and  the  link  throughput  and  interference  constraint  at  the  MAC  layer.  As 
shown  in  Sec.  ??  and  Sec.  V-D,  this  relationship  varies  with  the  application  type  (for  example, 
whether  acknowledgement  is  needed  to  complete  a  successful  data  transmission)  and  the  use  of 
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Fig.  3.  Spectrum  opportunity  detection  at  the  secondary  tx 


handshaking  signaling.  In  other  words,  it  depends  on  whether  the  roles  of  the  transmitter  and 
receiver  need  to  be  reversed  during  the  process  of  communicating  a  data  packet. 

IV.  Spectrum  Opportunity  Detection:  Figures  of  Merit 
In  this  section,  we  specify  the  figures  of  merit  for  spectrum  opportunity  detection. 

PHY  Performance.  The  spectrum  opportunity  detection  can  be  formulated  as  a  binary 
hypothesis  testing  problem.  Let  I(A,  d,  rx)  denote  the  presence  of  primary  receivers  within 
distance  d  to  the  secondary  transmitter  A,  and  I(B,  d.  tx)  the  presence  of  primary  transmitters 
within  distance  d  to  the  secondary  receiver  B.  Let  !(•,•,■)  denote  the  complement  of  I(-,  ■). 

The  two  hypotheses  are  given  by 

Hq  :  opportunity  exists,  i.e.,  I(A,  77,  rx)  n  l(B,  RF  tx). 

?ii  ■  no  opportunity,  i.e.,  I(A,  77,  rx)  U  I(B,  Rp  tx). 

The  figures  of  merit  at  the  physical  layer  are  given  by  the  probabilities  of  false  alarm  PF  and 
miss  detection  Pmd'- 

PF=Pr{ decide  Hx  \  H0},  PMD=Pr{ decide  H0  \  Hi}. 

The  performance  of  the  detector  is  specified  by  the  receiver  operating  characteristic  (ROC)  curve, 
which  gives  1  —  Pmd  (probability  of  detection  or  detection  power  denoted  by  PD)  as  a  function 
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of  PF.  See  Fig.  4  for  an  illustration.  In  general,  reducing  PF  comes  at  a  price  of  increasing 
Pmd  and  vice  versa.  The  tradeoff  between  false  alarm  and  miss  detection  is  thus  crucial,  and  the 
operating  characteristics  of  the  spectrum  sensor  should  be  designed  by  considering  the  impact  of 
detection  errors  on  the  MAC  performance.  As  a  consequence,  the  relationship  between  PHY  and 
MAC  needs  to  be  carefully  examined.  In  particular,  as  for  LBT  with  perfect  ears  (see  Fig.  3), 
the  tradeoff  between  false  alarm  and  miss  detection  can  be  achieved  by  adjusting  the  detection 
range  rD  e  (0 ,  ry  +  Rp\. 


Fig.  4.  ROC  curve  of  LBT  with  perfect  ears. 

MAC  Performance  The  MAC  layer  performance  is  measured  by  the  throughput  of  the 
secondary  user  and  the  interference  to  the  primary  users.  The  design  objective  is  to  maximize 
the  throughput  of  secondary  users  while  limiting  the  interference  to  the  primary  users. 

The  figures  of  merit  at  the  MAC  layer  are  thus  given  by  the  probability  P$  of  successful  data 
transmission  and  the  probability  Pc  of  colliding  with  primary  users1. 

Ps  =  Prjsuccessful  data  transmission},  (1) 

Pc  =  Pr{A  transmits  data  |  I(A,  ry,rx)}.  (2) 

'We  assume  that  collisions  with  primary  users  are  caused  by  data  transmissions.  We  ignore  the  interference  from  the 
transmission  of  acknowledgement  and  handshaking  signaling  such  as  RTS/CTS  due  to  their  short  duration. 
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Note  that  Pc  is  conditioned  on  I(A,  ry,rx)  instead  of  H\.  Clearly,  Pr[I(A,  77,  rx)]  <  Pr[7-C]. 
This  further  complicates  the  relationship  between  PMD  and  Pc  as  shown  in  Sec.  V. 

V.  Performance  of  LBT  with  Perfect  Ears 

In  this  section,  we  analyze  the  performance  of  LBT  with  perfect  ears  in  a  Poisson  distributed 
primary  network.  First,  since  false  alarm  and  miss  detection  are  conditioned  on  Hq  and  H\ 
respectively,  we  derive  the  expression  for  the  probability  of  spectrum  opportunity  Pr[Ho]. 
Second,  based  on  the  expression  for  Pr[Ho\,  we  obtain  the  analytical  expressions  for  PHY 
performance  measures  (false  alarm  probability  PF,  miss  detection  probability  Pmd )  and  MAC 
performance  measures  (successful  transmission  probability  P$,  collision  probability  Pc).  Last, 
we  compare  applications  requiring  guaranteed  delivery  with  those  relying  on  best  effort  (for 
example,  media  streaming  and  network  gaming),  and  also  study  the  impact  of  MAC  handshaking 
signaling  on  the  performance  of  LBT. 

A.  Probability  of  Spectrum  Opportunity  Pr[Hf\ 

Recall  the  definition  of  the  conditional  probability, 

Pr[Ho]  =  Pr{I(A,  r/,rx)  n  I(B,  f?/,tx)}, 

=  Pr{I(A,  r/,rx)  j  1(5,  f?/,  tx)}  •  Pr{II(5,  5/,  tx)}.  (3) 

Next  we  will  calculate  the  two  probabilities  in  the  above  expression  one  by  one. 

□  Pr{l(BJHtx)} 

Since  primary  transmitters  satisfy  a  Poisson  process  with  density  pX,  it  is  easy  to  see  that 

Pr{I(B,  Rh  tx)}  =  exp(—pAB),  (4) 

where  A B  =  Xn Rj. 


□  Pr{I(A,ry,rx )  j  I(B,  f?/,  tx)} 

Let  d  be  the  distance  between  A  and  B,  Sc(d,  rr,  r2)  denote  the  area  within  a  circle  with  radius 
?T  centered  at  A  but  outside  the  circle  with  radius  7*2  centered  at  B.  The  area  of  Sc(d,  rr,  rf) 
equals  7r r\  —  5/(d,  r±,  r2),  where  Sj(d,  ri,  r2)  is  the  common  area  of  two  circles  with  radius  ri 
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and  r2  and  centered  cl  apart  (see  Appendix  A  for  its  expression).  Then  by  considering  the  fact 
that  the  primary  receiver  within  77  of  A  can  only  communicate  with  the  primary  transmitter 
within  77  +  Rp  of  A  and  using  the  Total  Probability  Theorem,  we  have  that 

Pr{I(A,  77, rx)  | 

OO 

=  ^  Pr{I(A,  77, rx)  fl  ( k  tx  G  Sc(d,  rj  +  Rp ,  Ri))\  I( B ,  Rj,  tx)}, 

k=0 

OO 

=  ^  Pr{I(A,  77,rx)|  (k  tx  G  Sc(d,  77  +  Rp,  Rj))  D  I(B,  Rj,  tx)} 

k=0 

■  Pr{k  tx  G  Sc(d,  77  +  Rp,  Ri)}.  (5) 

Here  ”A:  tx  G  «Sc(d,  77  +  Rp,  Rj)”  means  that  there  exist  k  primary  transmitters  inside  the  region 
Sc(d,  77  +  Rp,Rj).  In  the  lase  step,  we  use  the  fact  that  the  spatial  distribution  of  primary 
transmitters  in  Sc(d,rj  +  Rp,  R.j)  is  independent  with  that  in  the  Rj  circle  of  B. 

Similarly  to  Pr{I(B,  Rj,  tx)},  Pr{k  tx  G  Sc(d,i'j  +  Rp,  f?/)}  is  given  by 

Pr{k  tx  G  Sc(d ,  77  +  Rp,  f?/)}  =  exp(-j7AAi),  (6) 

where  AAi  =  A Sc(d,  rj  +  Rp,  R /). 

Let  Pi  =  Pr{I(A,  r/,rx)|  (one  tx  G  Sc(d,rj  +  Rp,  R^)  fl  1(5,  A?/,  tx)},  then  since  each  pri¬ 
mary  transmitter  choose  its  receiver  independently,  it  follows  that 


Pr{I(A,  ?7,rx)|  (A;  tx  G  <Sc(d,  r/  +  Rp,  Rj))  fl  I(B,  Rj,  tx)} 

=  Pk={  ff  1 

J  J  $c(d,  rj  +  Rp,  A?/) 

yScCrf.rj+TJp,/?./) 

Here  r  is  the  radial  coordinate  in  the  polar  coordinate  system,  and  6  is  the  angular  coordinate. 
We  choose  the  secondary  transmitter  A  as  the  origin. 

□  Pr[H0 } 


1  - 


SjfaRp,  77) 


nm 


rdixW 


(7) 


By  combining  (7),  (6),  (5),  and  (4)  with  (3),  we  can  express  the  probability  Pr[7io]  of  hypothesis 
H0  in  the  following  form. 


Pr[H0]  =  exp[-p(AAi(l  -  Pi)  +  As)], 

/ 

=  exp  I  — pA 


Sr(r,  Rp,  ry) 


\ 


nRl 


rdrdO  +  nBj 


(8) 
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The  detailed  description  of  Sc(d,  rI-\-Rp,  Rj)  in  the  polar  coordinate  system  is  given  in  Appendix 
B.  Note  that  the  integrand  SlH^rH  docs  not  depend  on  the  angular  coordinate  9,  so  by 
integrating  with  respect  to  9  first,  we  can  reduce  the  double  integral  in  (8)  to  a  single  integral 
JqI+Rp  s'HrHriYr  ■  9(r )  dr,  where  (Hr)  is  a  function  of  r  and  determined  by  the  shape  of 
Sc(d,  rr  +  Rp ,  R[).  If  we  take  into  account  the  symmetry  of  Sc(d,  77  +  Rp,  Rf).  the  computation 
can  be  further  simplified. 

B.  PHY  Performance 

1)  False  Alarm  Probability  PF: 

For  LBT,  false  alarm  occurs  if  and  only  if  the  secondary  transmitter  A  detects  the  presence 
of  primary  transmitters  under  Ho  (see  Fig.  5).  In  this  case,  false  alarm  probability  PF  is  given 
by 

PF  =  Pr (I(A,  rDl  tx)  |  H0}, 

=  1  -  Pr{I(A,rD,tx)  |  H0}.  (9) 


Fig.  5.  An  illustration  of  false  alarm  (Under  Tto,  there  exist  primary  txs  within  td  of  A) 


Use  the  Bayer  Rule, 

Pr{I(A,rD,tx)  |  Ho} 


=  Pr{I(A,  rD,  tx)  |  I(A,r/,rx)  nl(5,i?7,tx)}, 

Pr{I{A,  n,  rx)  |  I(A,  rD,tx)  (T  I(B,  Rj,  tx)}  •  Pr{I(A,rD,\x)  fl  I(B,  f?/,  tx)} 

Pr[Ho\ 


(10) 
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Next  we  will  calculate  the  two  probabilities  in  the  numerator  of  the  above  expression  one  by 
one. 

□  Pr{I(A,  77>,tx)  fl  1(5,  5/,tx)} 


Similarly  to  Pr{I(B,  5/,  tx)},  Pr{l(A,  rD,  tx)  D  1(5,  72/,  tx)}  is  given  by 

Pr {I(A  rDl  tx)  n  1(5,  72/,  tx)}  =  exp(-pAAB),  (11) 

where  A Ab  =  A[7t (r2D  +  Rj)  -  5/(7,  rD,  5/)]. 

□  5r{I(A,r/,rx)  |  I(A,  775,  tx)  D  1(5,  RI:  tx)} 


Use  the  Total  Probability  Theorem  again, 

Pr{l(A,  r/,rx)  |  I(A,  775,  tx)  D  1(5,  5/,  tx)} 

OO 

=  PrfflA,  rj,  rx)  n  [7  tx  G  <SC(7, 77  +  5P,  5/)  —  <SA2]  |  I(A,  m,  tx)  fl  1(5,  5/,  tx)}, 

k= o 

OO 

=  ]T5r{I(Ar/,rx)  j  [7  tx  G  <SC(7,  r/  +  5P,  5/)  —  5A2]  H  I(A,  rD,  tx)  fl  1(5,  5/,  tx)} 
k= 0 

•  Pr{k  tx  G  <SC(7, 77  +  5P,  5/)  -  £42},  (12) 

where  5A2  =  Sr{d.  ry>  Rj)  fl5c(7,  ?7  +  5p,  5/).  The  expression  for  the  region  SA2  is  given 
in  Appendix  B.  In  the  last  step,  we  use  the  fact  that  the  spatial  distribution  of  primary 
transmitters  in  Sc(d,  rj  +  Rp.  5/)  —  SA2  is  independent  with  that  in  the  rD  or  77  +  Rp  circle 
of  A  and  the  Rj  circle  of  5. 

For  Pr{k  tx  G  Sc(d,  ri  +  Rp,  Rj)  —  SA2},  it  is  easy  to  see  that 

Pr{k  tx  G  Sc(d,  r  /  +  5P,  5/)  -  SA2}  =  — A'41  -  -  ^11  exp(-p(XAi  -  AA2)),  (13) 

where  AA2  =  A5A2. 

Let  52  =  5r{I(A, 77, rx)  |  [one  tx  G  Sc(d,rI  +  Rp,RI)  —  SA2]r\I(A,rD,tx)r\I(B,RI,tx)}, 

then  we  have  that 


□  PF 


Pr{I(A,  r/,  rx)  |  [7  tx  G  5C(7, 77  +  Rp,  Rj)  —  5A2]  n  I(A,  rD,  tx)  D  1(5,  5/,  tx)} 


( 

\ 

[[  1 

’  Sr(r,Rv,rry 

rdrcW 

1  JJ  5C(7,  r/  +  Rp,  Ri)  —  SA2 

nR2p 

\Sc(d,r i+Rp  ,Ri)—Sa2 

/ 
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Substitute  (13),  (14)  into  (12),  and  use  the  formula  exp(x)  =  Y1T=o  fr,  and  then  substitute 
(8),  (11)  and  (12)  into  (10), 

Pf  =  l- exp[-p((\A1- \2A2)(1- P2)  +  \Ab  ~  XAi(l- Pi)  -  \b)], 

S^r,  Rp,  77) 


=  1  —  exp 


-pX  7 tPd  -  Si(d,  I'D,  5/) 


Sa2 


-rdrd# 


•  (15) 


2 )  Miss  Detection  Probability  Pmd •' 

For  LBT  miss  detection  occurs  if  and  only  if  the  secondary  transmitter  A  does  not  detect  the 
presence  of  primary  transmitters  under  5o  (see  Fig.  6).  We  thus  have 

Pmd  =  Pr (I(A  rD,  tx)  |  PR}, 


=  Pr{I(A,  rD,  tx)  |  I(2l,77,rx)  U  1(5,  5/,tx)}. 


(16) 


Fig.  6.  An  illustration  of  miss  detection  (Under  Hi,  there  exists  no  primary  tx  within  rn  of  d) 


In  order  to  use  our  known  results  about  false  alarm  probability  PF,  we  write  Pmd  as  follows: 


P, 


MD  — 


Pr{I(A,  rDl  tx)}  —  Pr{I(A,  rD,  tx)  fl  I(A,  77,  rx)  D  1(5,  Rj.  tx)} 

5r[5t]  : 


Pr{I(A,  rDl  tx)}  —  Pr{I(A,  rD,  tx)  fl  I(A,  77,  rx)  fl  1(5,  5/,  tx)} 

1  -  Pr[H0]  ' 


(17) 


Since 


Pr{I(A,rD,tx)}  =  exp(-pAA3),  (18) 

Pr{I(A,  rD,  tx)  fl  I(A,  77,  rx)  fl  1(5,  5/,  tx)} 

=  Pr{I(A,  77,  rx)  |  1(^4,  775,  tx)  fl  1(5,  5/,  tx)}  •  Pr{I(A,  77?,  tx)  fl  1(5,  5/,  tx)},  (19) 
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where  A  .43  =  Xurf,  it  follow  from  (8),  (11),  and  (12)  that  the  expression  for  PMD  is  obtained 
as  below. 


Paid 


exp(—  pXTrrfj)  —  exp 

1  1 

1 

(nirjj+Rj)  +  ff  l(^’  l} rdrdO  J 

\  Sc(d,ri+Rp,Ri)-SA2  p  J 

1  —  exp 

-pA(  ff  SA^nlrdrd6  +  wRX 

{Sddpj+Rp’RP  p  ) 

{20) 


3)  Properties: 

□  Property  1  For  fixed  p\,  77,  d,  Rf,  and  Rp,  PF  ]  and  PMD  J.  as  rD  j. 


Proof.  Since 


7r r2D  -  Sr(d,  rD,  R /) 


<SA2 


SjjpRpP^ 
7T  Rp 


rdrdO 


<SA2 


Sr(  r\  Rp,  rj) 


nRl 


rdrdO 


and  the  area  SA2  of  the  double  integral  increases  as  rD  increases,  it  follows  from  the 
monotonicity  of  exponential  function  that  PF  is  an  increasing  function  of  rp. 

On  the  other  hand,  since 


Pmd 


exp(—  pXnrj-,)  / 1  —  exp 

pxInR2;  5/(d,rD,f?/)  +  ff  l]  rdrdO J 

\  Sc(d,rI+Rp,RI)-SA2  P  J 

1 

1  —  exp 

—pX  \  ff  <S/(r’p2,r/Vdrdfl  +  tt ) 

fs^n+RpPP  v  J  _ 

by  considering  the  monotonicity  of  exp(—  p\nr2D),  S[(d,  rD.  Rf),  and  SA2  with  respect  to 
rD,  we  conclude  that  PMD  is  a  decreasing  function  of  rp-  □ 


□  Property  2  For  fixed  rp,  rp  d,  Rp  and  Rp,  PF  |  as  pX  |. 

Proof.  Since  from  the  definition  of  SA2,  we  know  that  tt r2D  is  always  larger  than  or  equal 

to  Si(d,  rD,  f?/)  +  ff  Sl(r^fri\drd6 ,  based  on  the  monotonicity  of  exponential  function, 

SA2  ^  P 

we  can  conclude  that  PF  |  as  pX  j  given  other  parameters.  □ 


□  Property  3  For  fixed  rp,  rp  d,  Rp  and  Rp,  PMp  j  as  pX 
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Proof.  Let 


u  =  pX 
si  =  nr2D 

s2  =  7T R)  -  Sjid,  rD,  i?/)  + 


Si(r,  Rp,  r/) 


Sc(d,ri+Rp,Ri ) —Sa2 


ttR2 


rdrdd 


s  3  = 


Si(r ,  Rp,  r/) 


^Rl 


rdrd0  +  7ri?j 


Sc(d,ri+Rp,Ri) 

and  by  recalling  the  definition  of  Sc(d,  rj  +  R,p,  Rj  )  and  Sa2,  we  can  easily  find  that  u, 
si,  s-2,  and  s3  satisfy  the  following  properties:  all  of  them  are  positive,  and  si  +  s2  >  s3, 
s2  <  s3.  Then  we  only  need  to  show  the  function 

exp(— usi)(l  -  exp(— «s2)) 

JW  = - } - 7 - 7 - 

1  —  exp(— us3) 

is  a  decreasing  function  of  u  (u  >  0)  with  si,  s2,  s3  >  0  and  s\  +  s2  >  s3,  s2  <  s32. 

By  replacing  exp(— u)  by  x  (0  <  x  <  1),  we  obtain  a  new  function  of  x. 


rpS  1  (  1  /y»^2  ^ 

g(x)  =  — - ^  =  xSl~(S3~S2) 


x 


-S3  S2  _  53 


1  -  Xs3  1  -  Xs 3 

Since  x  =  exp(-u)  is  monotonically  decreasing  and  xSl_^3_S2)  (si  —  (s3  —  s2)  >  0)  is 
monotonically  increasing,  it  is  equivalent  to  show  that  h{x)  =  xS31^~:rS3  is  an  increasing 
function  of  x  for  0  <  x  <  1. 

Take  the  derivative  of  h(x). 


h'(x)  = 


x 


S3-S2-1 


[s3(l  —  Xs 2)  —  S2(l  —  Xs3) 


s  2 


-S3 


(1  -  xS3)2 

Next  we  will  show  that  s3(l  —  xS2)  —  s2(l  —  xS3)  >  0,  i.e.,  1~^2  >  1-^S3  for  0  <  s2  <  S3. 

To  prove  this  inequality,  we  define  a  function  of  a 

p(a )  =  - — (0<x<  1). 
a 

Take  the  derivative  with  respect  to  a, 

^  —  (ln(x))xaa  —  (1  —  xa) 

p  [a)  =  - 


“If  S2  =  S3,  then  f(u)  =  exp(— usi)  and  the  proof  is  trivial 
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Denote  the  numerator  of  p'(a)  by  q(a),  then  we  have  that 

q( 0)  =  0,  q\a)  =  — ( ln(x))2xaa  <  0,  for  all  a  >  0 
=>•  q(a)  <  0,  for  all  a  >  0,  i.e.,  p'{a)  <  0,  for  all  a  >  0 

So  p(a)  is  monotonically  decreasing  with  respect  to  a.  It  follows  that  the  inequality  holds 
for  0  <  s2  <  s3.  Hence,  h(x)  is  an  increasing  function  of  x  (0  <  x  <  1).  □ 


C.  MAC  Performance 

1)  Collision  Probability  Pc: 

Since  the  secondary  transmitter  A  transmits  data  if  and  only  if  A  detects  no  nearby  primary 
transmitters  (see  Fig.  7),  by  recalling  (1),  we  have 

Pc  =  Pr{I(A,rD,tx)  |  I(A,r/,rx)}.  (21) 

Note  that  collision  probability  Pc  does  not  depend  on  the  activities  of  primary  users  near  the 
secondary  receiver  B.  Moreover,  By  comparing  (21)  with  (16),  we  can  easily  see  that  Pc  f  Pmd- 


Fig.  7.  An  illustration  of  collision  (there  exists  no  tx  within  ro  of  A  =£•  Collision  with  some  rx  within  n  of  A) 


Based  on  the  basic  property  of  conditional  probability,  we  have  that 

_  Pr{I(A,rDltx)  fl  1(^4, rj,rx)} 

C  Pr{I(A,  rj,rx)} 


(22) 
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Since  according  to  the  Displacement  Theorem  primary  receivers  admits  a  spatial  Poisson 
distribution  with  density  pX,  it  follows  that 

Pr{T(A,  ry,rx)}  =  1  —  exp(—  pXrcr]).  (23) 

Then  by  using  the  similar  techniques  in  the  derivation  of  the  expression  for  Pr[Ho],  we  can 
obtain  the  expression  for  the  probability  in  the  numerator  of  (22).  So  collision  probability  Pc 
is  given  by 

exp(—pXnr2D)[l  -  ex-p(-pXn(r2T  -  I(rn,rJ,Rv)))\ 

f c  —  - : - 7 - 7 — 2\ - > 

1  —  exp(—  pXnrf 

where  I(rD,  ry,  Rp)  =  J^D  2r‘S/(^Rp)  dr.  If  we  substitute  the  expression  for  Sj(r,  ri,Rp)  into 
the  above  integral,  we  can  actually  obtain  a  closed-form  expression  for  /(?’D,ry,  Rp)  (see  Ap¬ 
pendix  C). 

2 )  Successful  Transmission  Probability  P$: 

For  applications  with  guaranteed  delivery,  an  acknowledgement  (ACK)  signal  from  B  to  the 
secondary  transmitter  A  is  required  to  complete  a  data  transmission.  In  this  case  successful 
data  transmission  occurs  if  and  only  if  the  ACK  signal  is  successfully  received  at  A.  But  since 
spectrum  opportunity  is  asymmetric  (see  Sec.  Ill),  correctly  identified  opportunity  may  not  lead 
to  a  successful  transmission  (see  Fig.  8).  On  the  other  hand,  as  shown  in  Fig.  9,  if  there  are 
some  primary  receivers  within  ry  of  A,  but  no  primary  transmitters  within  If  of  A  and  If  then 
miss  detection  leads  to  a  successful  transmission. 


Fig.  8.  Successful  data  transmission  but  failed  ACK 
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S  S  Primary  Tx 
•  Primary  Rx 


Fig.  9.  Miss  detection  leads  to  successful  transmission. 


Let  the  effective  range  rE  =  max{rD,  Rj},  then  followed  from  the  above  facts,  we  can  express 
the  probability  Ps  of  successful  data  transmission  as  follows: 

Ps  =  Prjsuccessful  data  transmission}, 

=  Pr{successful  transmission  |  H0}  ■  Pr[H0]  +  Prjsuccessful  transmission  |  Hi}  ■  Pr[H\ ], 
=  Pr[I(A,  rEl  tx)  |  Ho]  ■  Pr[H0]  +  Pr[I(A,  rE,  tx)  n  I(B,  Rj,  tx)  |  Hi]  ■  Pr[Hi].  (25) 

By  comparing  (25)  with  (9)  and  (16),  we  can  easily  see  that 

Pr{successful  transmission  |  Ho}  <  1  —  Pfi 
Pr{successful  transmission  |  Hi}  <  Paid- 

By  using  the  definition  of  Ho  and  Hi,  we  can  simplify  the  expression  for  Ps  into  the  following 
form: 

Ps  =  Pr{I(A,rEl  tx)  D  I(P,  P/,  tx)}, 

=  exp[— pA(7r(r|  +  Rj)  -  5/(d,  rE,  Ri ))]•  (26) 

From  (26),  it  seems  that  the  success  probability  Ps  does  not  depend  on  77.  But  since  the  collision 
probability  Pc  (see  (24))  depends  on  77,  and  there  is  a  constraint  on  Pc  ( Pc  <  C)>  it  follows 
that  given  different  77,  rD  probably  needs  to  be  chosen  differently  to  meet  the  constraint  on  Pc 
and  therefore  it  results  in  different  P$. 

We  now  consider  the  success  probability  Ps  for  applications  relying  on  best-effort  delivery 
(for  example,  media  streaming  and  network  gaming).  For  these  applications,  ACKs  are  not 
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necessary.  Since  ACK  is  not  involved  in  opportunity  detection  and  collision  with  primary  users, 
given  the  same  parameters  (p,  A,  Rp,  Rj,  d,  ry,  and  rD)  we  have  the  same  expressions  for  false 
alarm  probability  PF,  miss  detection  probability  Pmd,  and  collision  probability  Pc  as  before 
but  different  successful  transmission  probability  Ps. 

For  best-effort  delivery,  data  transmission  is  successful  if  and  only  if  data  is  received  success¬ 
fully  at  the  secondary  receiver  B.  So  the  probability  Ps  of  successful  data  transmission  is  given 
by 

Ps  =  Pr {successful  data  transmission}, 

=  Pr{successful  transmission  |  7 Y0}  •  Pr[H0]  +  Pr{successful  transmission  |  Hi}  •  Pr[Hi\, 
=  Pr[I(A,rD,tx)  |  H0]  ■  Pr[H0]  +  Pr[I(A,  rD,  tx)  n  I(B,  RI:  tx)  |  Hf\  ■  Pr[Hi], 

=  Pr{I(A,  r_o,  tx)  fl  I(B,  Rj,  tx)}, 

=  exp[— pA(7r(rfj  +  Pj)  —  Si(d,rD,  Ri))\.  (27) 

By  comparing  (27)  with  (9)  and  (16),  we  can  easily  see  that  for  best-effort  delivery: 

Pr{successful  transmission  j  Hq}  =  1  —  PF, 

Pr{successful  transmission  |  Hi}  <  Pmd- 

D.  Performance  of  RTS/CTS  enhanced  LBT 

The  fundamental  deficiency  of  LBT  resembles  the  hidden  and  exposed  terminal  problem  in 
the  conventional  ad  hoc  networks  of  peer  users.  It  is  thus  natural  to  consider  the  use  of  RTS/CTS 
handshaking  signaling  to  enhance  the  detection  performance  of  LBT. 

For  RTS/CTS  enhanced  LBT,  spectrum  opportunity  detection  is  performed  jointly  by  the 
secondary  transmitter  A  and  the  secondary  receiver  B  through  the  exchange  of  RTS/CTS  signals. 
The  detailed  steps  of  RTS/CTS  enhanced  LBT  is  listed  as  below. 

•  A  detects  primary  transmitters  within  distance  td-  If  it  detects  none,  A  sends  B  a  Ready- 
to-Send  (RTS)  signal. 

•  If  B  receives  the  RTS  signal  from  A  successfully,  which  automatically  indicates  the  absence 
of  primary  transmitters  within  distance  R  j,  then  B  replies  with  a  Clear-to-Send  (CTS)  signal. 

•  Upon  receiving  the  CTS  signal,  A  transmits  data  to  B. 
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Then  we  will  derive  the  expressions  for  false  alarm  probability  PF,  miss  detection  probability 
Pmd,  collision  probability  Pc,  and  successful  transmission  probability  Ps  one  by  one. 

1)  False  Alarm  Probability  PF: 

For  RTS/CTS  enhanced  LBT,  the  observation  space  of  opportunity  detection  comprises  the 
RTS  and  CTS  signals.  So  false  alarm  occurs  if  and  only  if  the  RTS/CTS  exchange  fails  under 
Ho  (see  Fig.  10  and  Fig.  11).  It  follows  that  false  alarm  probability  PF  is  given  by 

PF  =  Pr{I(A,rE,tx)  |  H0}  (28) 


By  comparing  (28)  with  (9),  we  find  that  they  are  almost  the  same  except  that  rD  is  replaced 
by  rE  in  (28).  Hence  we  have  that 


PF  —  1  —  exp 


— pX 


Si(d,  rE,  Ri) 


S/(r,  f^yr/) 


TtRp 


rdrd# 


(29) 


where  SA \  =  S,id.  rE,  fZ/)  D  Sr(d.  ry  +  Rp.  Rf  ),  and  the  detailed  form  of  SA ]  can  be  found  in 
Appendix  B. 


Fig.  10.  Case  1:  A  detects  primary  transmitter. 

2 )  Miss  Detection  Probability  Pmd- 

For  RTS/CTS  enhanced  LBT,  miss  detection  occurs  if  and  only  if  the  RTS/CTS  exchange  is 
successful  under  Hi  (see  Fig.  12).  So  miss  detection  probability  Pmd  is  given  by 


Paid  =  Pr{I(A,  rD,  tx)  fl  I (B,  RI}  tx)  D  I(A,  Rj,  tx)  |  Hi}. 


(30) 
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(b)  there  exists  tx  within  Rj  of  A 


Fig.  11.  Case  2:  A  fails  to  receive  CTS. 


(a)  A  detects  no  primary  tx  (b)  B  receives  RTS 

=>■  no  tx  <  td  of  A  =>-  no  tx  <  Rj  of  B 


(c)  A  receives  CTS 
=>■  no  tx  <  Rj  of  A 


Fig.  12.  An  illustration  of  miss  detection 


Then  by  using  Bayes  rule,  we  can  express  (30)  in  the  form  of  some  know  probabilities. 


Pmd  = 


Pr{I(A,  te,  tx)  Dl(v4,rj,rx)  nI(B,i?/,tx)} 
Pr\n  i]  ! 


^1  —  Pr{I(T,r/,rx)  |  I(A,  r£,tx)  nl(f?,f?/,tx)}^  ■  Pr{I(A,  rs,tx)  n  I(f?,  Ri,  tx)} 


1  -  Pr[H0 


(31) 
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By  recalling  (8),  (11),  and  (12),  we  have  that 

=  exp(—  p\AB)[l  -  exp(— p(A^ti  -  A\2)(l  -  P2))} 
MD  1  -  exp[-p(\A1(l  -  P±)  +  \B)} 


exp[— pA(7r(r|  +  i?.f)  -  Si(d,  rE,Ri ))] 

1  —  exp  |  —  pX  JJ  Si(rjip,ri)  rdrdg  j 

V  Sc(d,n+Rp,Ri)-SA  i  P  Jl 

1  —  exp 

-pX(  u 

\Sc(d.rj+i?p.Rr)  P  ) 

where  A^2  =  ASai- 
3)  Collision  Probability  Pc,: 

For  RTS/CTS  enhanced  LBT,  given  that  there  exists  some  primary  receiver  within  the  inter¬ 
ference  range  ry  of  the  secondary  transmitter  A,  collision  occurs  if  and  only  if  the  RTS/CTS 
exchange  is  successful  (see  Fig.  13).  So  collision  probability  Pc  is  given  by 

Pc  =  Pr{I(A,  rE,  tx)  n  I(B,  RT,  tx)  |  I(A4,  ry,  rx)}.  (33) 


(a)  A  detects  no  primary  transmitter 
sexists  no  tx  within  m  of  A 


and 


(c)  A  receives  CTS  successfully 
sexists  no  tx  within  Ri  of  A 


(b)  B  receives  RTS  successfully 
sexists  no  tx  within  Rj  of  B 


(d)  A  sends  data  and  causes  interference 


Fig.  13.  An  illustration  of  collision 
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By  comparing  (33)  with  (16),  we  have  that 


Pc  = 


Pr[Hi 


-Pmd, 


Pr{I(A,  77,  rx)} 
exp[-pA(7r(r|  +  Rj)  -  Si{d,  rE,  Ri ))] 


1-«P  -pA  ff 

\  Sc(d,ri+Rp,Ri)—SAi 


rdrd^J 


1  —  exp  ( pXnrj ) 


(34) 


Note  that  Pr[Hf\  =  Pr{I(A,rj,rx)  U  I(B,  i?/,  tx)}  >  Pr{I(A,  77,  rx)},  so  Pc  >  Pmd- 
4)  Successful  Transmission  Probability  P$: 

For  RTS/CTS  enhanced  LBT,  a  data  transmission  is  successful  if  and  only  if  the  RTS/CTS 
exchange  is  successful  (see  Fig.  14).  The  secondary  transmitter  A  sends  data  to  the  secondary 
receiver  B  only  if  A  receives  CTS  from  B.  Since  before  B  replies  with  CTS  B  must  receive 
RTS  successfully,  B  can  also  receive  data  successfully.  So  the  probability  Ps  of  successful  data 
transmission  is  given  by 


Ps  =  Pr{I(A,rE,tx)  nl(S,i?/;tx)}, 

=  exp[— pA(7r(r|  +  R ])  -  Sj{d,  rE,  Ri ))].  (35) 


Notice  that  Ps  of  RTS/CTS  enhanced  LBT  is  identical  to  that  of  LBT  without  RTS/CTS  for 
guaranteed  delivery  in  (27).  Moreover,  using  RTS/CTS  enhanced  LBT,  the  throughput  is  the 
same  for  guaranteed  delivery  and  best-effort  delivery. 


(a)  A  detects  no  primary  transmitter 
=>exists  no  tx  within  rE  of  A 


(b)  B  receives  RTS  successfully  (c)  A  receives  CTS  successfully 
Sexists  no  tx  within  Rj  of  B  Sexists  no  tx  within  II 1  of  A 


Fig.  14.  An  illustration  of  successful  transmission 


5 )  Properties:  By  using  techniques  similar  to  proving  these  properties  for  LBT,  we  can  show 
that  Property  1,2, 3, 4  also  hold  for  RTS/CTS  enhanced  LBT. 
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VI.  Numerical  Examples 

In  this  section,  we  first  study  the  PHY  layer  performance  of  opportunity  detection  for  LBT  by 
observing  the  ROC  curves,  and  then  we  investigate  the  impact  of  application  and  MAC  hand¬ 
shaking  signaling  on  the  translation  from  PHY  layer  performance  to  MAC  layer  performance. 
All  the  figures  in  this  section  are  plotted  based  on  the  expressions  derived  in  Sec.  V. 

A.  PHY  Performance 


Q 

0. 


p 


F 


Fig.  15.  ROC  performance  comparison  (p  =  0.03,  A  =  10/2002,  Rp  =  200 m,  Ri  =  250 m,  d  =  200 m,  ri  =  200m/0.9). 


In  Fig.  15,  we  plot  the  ROC  curves  for  both  LBT  and  RTS/CTS  enhanced  LBT,  and  we  can 
observe  the  following  facts: 

•  (0,0)  does  not  belong  to  the  ROC  curve  of  RTS/CTS  enhanced  LBT.  This  is  due  to  the 
fact  that  the  effective  detection  range  is  bounded  above  Ri,  since  to  receive  the  CTS  signal 
successfully,  there  cannot  be  primary  transmitters  within  Rj  of  A.  In  other  words,  a  detection 
range  rD  <  If  leads  to  the  same  (If,  Pmd )  as  m  =  Ri- 

•  The  ROC  performance  of  RTS/CTS  enhanced  LBT  is  always  better  than  or  equal  to  that  of 
LBT  when  rD  >  Rf.  When  rD  <  Rj,  we  can  use  randomized  detection  with  parameter  pT 
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to  complete  the  irregular  ROC  curve  of  RTS/CTS  enhanced  LBT  in  this  range.  Specifically, 
if  the  RTS/CTS  exchange  succeeds,  we  decide  H0;  otherwise,  if  the  RTS/CTS  exchange 
fails,  we  decide  Ho  with  probability  pT  and  H\  with  probability  1  —  pT.  In  this  case,  there 
is  a  straight  line  that  connects  (0,0)  with  the  point  (PF,PD)  where  rD  =  Rj,  i.e.,  the 
starting  point  in  Fig.  15.  It  is  easy  to  see  that  in  this  range,  the  ROC  performance  of  this 
randomized  LBT  with  RTS/CTS  may  be  worse  than  that  of  LBT. 


B.  Impact  of  Application 

Since  the  design  objective  is  to  maximize  Ps  (representing  link  throughput)  under  a  collision 
constraint  Pc  <  C>  we  show  Ps  as  a  function  of  Pc  in  Fig.  16. 


Fig.  16.  Success  probability  vs  interference  constraint  (p  =  0.03,  A  =  10/2002,  Rp  =  200 m,  Ri  =  250 m,  d  =  200m, 
n  =  200m/0.9). 


From  Fig.  16,  we  observe  that  even  though  the  detection  performance  at  the  physical  layer 
is  the  same  (see  Fig.  15),  the  MAC  layer  performance  can  be  different  depending  on  the 
applications.  When  the  collision  constraint  is  tight,  the  throughput  is  the  same  for  these  two  types 
of  applications.  The  collision  constraint  (  has  a  critical  value  Co  above  which  the  throughput  for 
best-effort  delivery  is  higher  than  that  for  guaranteed  delivery.  Fig.  17  shows  Co  as  a  function 
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of  the  primary  traffic  load  pX  (or  the  density  of  active  primary  transmitters).  We  can  see  that 
Co  is  a  decreasing  function  of  pX.  This  suggests  that  primary  systems  with  heavy  traffic  is  more 
suitable  for  spectrum  overlay  with  best-effort  delivery  applications. 


Fig.  17.  Critical  value  of  £  vs.  primary  traffic  load  (A  =  10/2002,  Rp  =  200 m,  m  =  Ri  =  250 to,  d  =  200m,  ri  = 
200?n/0.9). 

C.  Impact  of  MAC  Handshaking  Signaling 

Fig.  18  shows  Ps  as  a  function  of  Pc  for  LBT  in  two  typical  applications  (guaranteed  delivery, 
best-effort  delivery)  and  RTS/CTS  enhanced  LBT.  Note  that  using  RTS/CTS  enhanced  LBT,  the 
throughput  is  the  same  for  guaranteed  delivery  and  best-effort  delivery.  We  can  see  that  at  the 
RTS/CTS  enhanced  LBT  may  lead  to  lower  throughput  when  the  collision  constraint  is  loose 
and  the  application  relies  on  best-effort  delivery.  Note  that  even  random  detection  is  used  in 
RTS/CTS  enhanced  LBT  for  rp  <  Rj  to  complete  the  ROC  curve,  its  throughput  can  still  be 
lower  than  that  of  LBT  for  best-effort  delivery  when  (  is  relatively  large  as  shown  in  Fig.  19. 

VII.  Conclusion 

In  this  report,  by  considering  a  Poisson  distributed  primary  network,  we  have  derived  the 
analytical  expressions  for  opportunity  detection  performance  at  both  PHY  and  MAC  layers. 
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Fig.  18.  Throughput  comparison  (p  =  0.03,  A  =  10/2002,  Rp  =  200 m,  Ri  =  250m,  d  =  200 m,  n  =  200m/0.9). 


c 


Fig.  19.  Throughput  comparison  ( p  =  0.03,  A  =  10/2002,  Rp  =  200m,  Ri  =  250m,  d  =  200 m,  n  =  200?n/0.9). 
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Based  on  these  expressions,  we  illustrate  the  complex  dependency  of  the  relationship  between 
PHY  and  MAC  on  the  application  type  (guaranteed  delivery  vs.  best-effort  delivery)  and  the 
use  of  handshaking  signaling  such  as  RTS/CTS  at  the  MAC  layer  by  some  numerical  examples. 
Specifically,  when  the  interference  constraint  is  tight,  we  should  adopt  RTS/CTS  handshaking 
to  improve  the  link  throughput,  but  RTS/CTS  handshaking  leads  to  decreased  throughput  when 
the  interference  constraint  is  loose.  Moreover,  in  primary  networks  with  relatively  heavy  traffic 
best-effort  delivery  applications  are  more  suitable  than  guaranteed  delivery  applications. 

Appendix 

A.  Intersecting  Area  Sj  of  Two  Circles 

Assume  that  the  radii  of  two  circles  are  r\  and  r2  respectively  and  the  distance  between 
the  centers  of  two  circles  is  d.  Without  loss  of  generality,  we  also  assume  that  r2  >  r\.  The 
expression  for  Sj  depends  on  the  relation  between  d  and  r2  —  r, ,  r1  +  r2  (see  Fig.  20),  and  it 
is  given  in  the  following  three  cases. 


Fig.  20.  Intersecting  area  of  two  circles 


□  Case  1  (0  <  d  <  r2  —  ry): 


Si(d,ri,r2)  =  Ttr\ 
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□  Case  2  (r2  —  rv  <  d  <  r\  +  r2)  [4]: 


Si(d,r1,r2) 


□  Case  3  (d  >  +  r2): 


r0  cos 


-i 


+  d2 


2  d  r2 


+  ?’2  cos  1 


+  d 2 


2  d  ri 


2 


S/(d,ri,r2)  =  0 


/f.  Sc(d,7'i  +  Rp,  Ri),  SA i  and  Sa2  in  the  Polar  Coordinate  System 
Here  we  pick  the  secondary  transmitter  A  as  the  origin. 

1 )  Sr(d.  ri  +  Rp,  Ri)  -  the  area  within  a  circle  with  radius  ri  +  Rp  centered  at  A  but  outside 
the  circle  with  radius  Ri  centered  at  B: 

>  Case  1:  Ri  >  d  (see  Fig.  21) 


Sc(d,  ri  +  Rp,  Ri) 


=  < 


Ri  —  d  <  r  <  rj  +  Rp,  9o(r)  <  6(r)  <27 r  —  9o(r), 


Ri  —  d  <  r  <  rj  +  Rr, 


9o(r)  <  9{r)  <2tt  —  90(r),  for  r  <  Ri  +  d 
0  <  9(r)  <  27 r,  otherwise 


where  6*o(r)  =  arccos  1  ^4 — 


2  dr 


>  Case  2:  Rj  <  d 


if  n  +  Rp  <  Rj  —  d, 

if  Ri  —  d  <  ri  +  Rp  <  Ri  +  d, 


if  Vi  -f-  Rp  <  Ri  -t-  d, 


Sc(d,  ri  +  Rp,  Ri) 


< 


0  <  r  <  ri  +  Rp,0  <  9  <  2tt 

!0  <  9(r)  <  2n, 

9o(r)  <  9(r)  <2i r  -  90(r), 
9o(r)  <  9{r)  <  2tt  -  90(r), 
0  <  9(r)  <  2tt, 


for  r  <  d  —  Ri 
otherwise 

for  d  —  Rf  <  r  <  Rr  +  d 
otherwise 


if  n  +  Rp  <  d  —  Ri, 

if  d  —  Ri  <  n  +  Rp  <  Ri  +  d, 


if  ri  +  Rp  >  Ri  +  d. 
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Fig.  21.  Illustration  of  Sc  (Case  1) 


2)  SA 1  =  Sc(d,  rE,  Rj)  fl  Sc(d,  rT  +  Rp,  R^: 

>  Case  1:  i?/  >  d,  d  +  R,j  <  77  +  Rp  (see  Fig.  22) 


Sai  —  < 


Ri  —  d  <  r  <  Ri,  9o(r)  <  9(r)  <27 r  —  9q (r) 

Ri  —  d  <r  <  fo,  0o (t)  <  9(r)  <2n~  9o(r) 

\  9o(r)  <  9(r)  <  2n  —  0o(t),  for  r  <  Ri  +  d 
Rj  —  d  <  r  <  re,  < 

I  0  <  9(r)  <  2ir,  otherwise 

Sc(d,rE,  Ri), 


if  rD  <  Ri, 

if  Ri  <  ro  <  Ri  +  d, 


if  rD  >  n  +  Rp- 


>  Case  2:  i?/  >  d,  d  +  R,j  >  77  +  Rp 

{Ri  —  d  <  r  <  Ri,  9q(t)  <  9(r)  <  2tt  —  90(r )  if  rE  <  Ri, 

Ri  —  d  <  r  <  rE,  90 (r)  <  9(r)  <2i r  —  90(r)  if  Ri  <  rD  <  n  +  Rp, 

Sc(d,  rE,  Ri),  if  rD  >  n  +  Rp. 
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Fig.  22.  Illustration  of  <Sai  (Case  1) 


>  Case  3:  |  <  Rj  <  d,  d  +  i?/  <  77  +  Rp 


—  < 


0  <  r  <  Rj, 

0  <r<  td, 


0  <  6(r)  <  27 r, 


for  r  <  d  —  Ri 


0o(r)  <  9{r)  <2ti  —  6o(r),  otherwise 


0  <  0(r)  <  2n, 


for  r  <  d  —  Ri 


0  <  0(r)  <  27r, 


0  <  r  <  rr>, 

Sc(d,rE,Ri), 

>  Case  4:  ^  <  Rj  <  d,  d  +  Rr  >  rr  +  Rp 

0  <  r  <  Ri, 


9o(r)  <  0(r)  <27 r  —  9o(r),  otherwise 

9o(r)  <  0(r)  <  2tt  —  90(r),  for  d  —  Ri<r<Ri+d 


otherwise 


if  r'n  <  Ri 

if  Ri  <  rE  <  Ri  +  d 

if  Ri  +  d  <  rD  <  n  +  Rp 
if  rD  >  n  +  Rp 


tSyll  —  < 


0  <  r  <  td, 

Sc(d,rE,  Ri), 


0  <  9{r)  <  27 r, 

for  r  <  d  —  Ri 

if  rD  <  Ri 

9o (r)  <  0(r)  <2n  -  90(r), 

otherwise 

0  <  9(r)  <  2n, 

for  r  <  d  —  Ri 

5 

if  Ri  <  rD  <  n  +  R1 

0o (0  <  0(r)  <  27T  -  0o(r-), 

otherwise 

if  rD  >  ri  +  Rp 
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t>  Case  5:  R.j  <  ^,  d  +  Rr  <  77  +  Rp 


Sai  =  { 


0  <  r  <  Ri,  0  <  9  <  2ir, 

0  <  r  <  I'd,  0  <  9  <  2i r, 

0  <  9(r)  <  27 r, 


0  <  r  <  I'd, 

0  <  r  <  I'd, 

Sc(d,rE,Ri), 


9q(t)  <  9(r)  <  2tt 
9o(r)  <  9(r)  <  27r 
0  <  9(r)  <  2i r, 


for  r  <  d  —  Ri 

5 

9o(r),  otherwise 
9o(r),  for  d  —  Ri  <  r  <  Ri  +  d 
otherwise 


if  I'd  <  Ri 

if  Ri  <  t_d  <  d  —  Ri 

if  d  —  Ri  <  rE  <  Ri  +  d 

if  Ri  +  d  <  I'd  <  ti  +  Rp 
if  I'd  >  ri  +  RP 


>  Case  6:  Ri  <  %,  d  +  Ri  >  77  +  Rp 

0  <  r  <  Ri,  0  <  9  <  2n, 

0  <  r  <  I'd,  0  <  9  <  2i r, 

<5,41  =  {  [  0  <  9(r)  <  2i r, 


for  r  <  d  —  Ri 
9q (r)  <  9(r)  <2n  —  9q{t),  otherwise 


3)  SA2  =  <Sc(d,  rD,  Ri)  n  Sc(d,  rr  +  Rp,  RT): 
[>  Case  1:  i?/  >  d,  d  +  Rj  <  77  +  Rp 


if  I'd  <  Ri 

if  Rf  <  I'd  <  d  —  Ri 

if  d  —  Ri  <  rD  <  ri  +  Rp 
if  I'd  >  ri  +  Rp 


<Sa2  —  \ 


Ri  —  d  <  r  <  I'd,  9o(r)  <  9(r)  <  2n  —  9o(r) 

%{r)  <  9(r)  <  2n  —  9o(r),  for  r  <  Ri  +  d 


Ri  —  d  <  r  <  I'd, 


Sc(d,rE,  Ri), 


0  <  9(r)  <  2n, 


otherwise 


if  I'd  <  Ri  —  d 

if  Ri  —  d  <  I'd  <  Ri  +  d 

if  Ri  +  d  <  I'd  <ri  +  Rp 
if  rD  >  1  'i  +  Rp 


>  Case  2:  i?/  >  d,  d  +  Rj  >  r/  +  Rp 

{4>  if  I'd  <  Ri  —  d 

Ri  —  d  <  r  <  rE,  90(r)  <  9(r)  <2n  —  9o(r)  if  Rj  —  d  <  rE  <  ri  +  Rp 
SAi,  if  rD>ri+Rp 


>  Case  3: 


Sa2  —  \ 


Ri  d  R  j  7  /  /  Rp 

0  <  r  <  ro,  0  <  9  <  2i r, 

!0  <  9(r)  <  2i r, 
90(r)  <  9(r)  <  2tt 
°o (r)  <  9{r)  <  2tt 
0  <  9(r)  <  2i r, 


y  Sc(d,rE,Ri), 


for  r  <  d  —  Ri 
9o(r),  otherwise 
9o(r),  for  d  —  Ri  <  r  <  Ri  +  d 
otherwise 


if  i'd  <  d  —  Ri 

if  d  —  Ri  <  rE  <  Ri  +  d 

if  Ri  +  d  <  I'd  <  ri  +  Rp 

if  i'd  >  ri  +  Rp 
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t>  Case  4:  f?/  <  d,  cl  +  Ri  >  r7  +  Rp 


$A2  —  \ 


0  <  r  <  rE, 
0  <  r  <  rE, 


0  <  9  <  2tt, 
j  0  <  9{r)  <  27 r, 

)  9o(r)  <  9(r)  <  27t 


[  Sc(d,rE,Ri), 


for  r  <  d  —  Ri 
9g(r),  otherwise 


if  rD  <  d  —  Ri 

if  d  —  Ri  <  rE  <  ri  +  Rp 

if  rD  >ri  +  Rp 


C.  Expression  for  l(rD.  rr.  Rp) 
Recall  that 


/(rD,r7,i?p)  = 


rrD«„Si(r>ri>Rv) 


2  r- 


dr 


By  using  the  expressions  of  Si(r,  ri,Rp)  in  Appendix  A  and  figuring  out  the  integral  in  the 
above  expression,  we  can  obtain  the  closed-form  expression  as  below. 

□  Case  1  (rj  <  Rv): 

•  rD  <  Rp-  n 


/(rD,r7,f?p)  = 


r2  2 
'I'D 

R2 

1  Lp 


rD  >  RP  +  r7 


I(tD,  rj,  Rp)  =  r] 


•  R  p  —  Tj  <  Tj)  <  Rp 


ri 


/(rD,r/,f?p 


±  2  :  '  D 

-r7  d - arccos 

2  7 r 


Rl  +  r2D  -  r2 


2„2 


rfr 


I'  D 


TT  Rp 


arccos 


2  Rpt  d 
rj  +  rD 


2  £  -2  Ri 


2  rirD 


- arcsin 

7T 


rf  +  R,i  -r2D\  r2rj  +  rf  +  /? 


2riRp 


□  Case  2  (r7  >  Rp): 


4:71  R2 


■yJ(rI  +  Rp  +  rD){rI  +  Rp-  rD)  (r/  -  f?p  +  rD)  -  r7  +  rD) 


rD  <  r7  -  f?p 


I(rD,rIlRp)  =  r2D 
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•  rD  >  n  +  Rp 

I(rD,  rI:  Rp)  =  r] 

•  rj  -  Rp  <  rD  <  n  +  Rp 

The  same  as  Case  1  (Rp  —  rj  <  ro  <  RP  +  rr) 
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